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Abstract. We have obtained high resolution, high signal- 
to-noise ratio CCD echelle spectra of 10 bright red giants 
in 3 globular clusters (47 Tuc, NGC 6752 and NGC 6397) 
roughly spanning the whole range of metallicities of the 
galactic globular cluster system. The analysis of this newly 
acquired material reveals no significant evidence of star- 
to-star variation of the [Fc/H] ratio in these three clusters. 
Moreover, a large set of high quality literature data (equiv- 
alent widths from high dispersion CCD spectra) was re- 
analyzed in an homogeneous and self-consistent way to in- 
tegrate our observations and derive new metal abundances 
for more than 160 bright red giants in 24 globular clusters 
(i.e. about 16% of the known population of galactic glob- 
ulars). This set was then used to define a new metallicity 
scale for globular clusters which is the result of high qual- 
ity, direct spectroscopic data, of new and updated model 
atmospheres from the grid of Kurucz (1992) , and of a care- 
ful fine abundance analysis; this last, in turn, is based on 
a common set of both atomic and atmospheric parameters 
for all the stars examined. Given the very high degree of 
internal homogeneity, our new scale supersedes the offsets 
and discrepancies existing in previous attempts to obtain 
a metallicity scale. The internal uncertainty in [Fe/H] is 
very small: 0.06 dex (24 clusters) on average, and can be 
interpreted also as the mean precision of the cluster rank- 
ing. Compared to our system, metallicities on the widely 
used Zinn and West's scale are about 0.10 dex higher for 
[Fe/H]> -1, 0.23 dex lower for -1 <[Fe/H]< -1.9 and 
0.11 dex too high for [Fe/H]< -1.9. The non-hnearity of 
the Zinn and West's scale is significant even at 3(7 level. A 
quadratic transformation is given to correct older values 
to the new scale in the range of our calibrating clusters 
(—2.24 <[Fe/H]ziv < —0.51). A minor disagreement is 
found at low metallicities between the metallicity scale 
based on field and cluster RR Lyrae variables {via a new 



calibration of the AS index) and our new cluster metallici- 
ties. It could be tentatively ascribed to non-linearity in the 
[Fe/H]— AS relationship. The impact of new metallicities 
on major astrophysical problems is exemplified through 
a simple exercise on the Oosterhoff effect in the classical 
pair M 3 and M 15. 

Key words: Stars: abundances - Stars: Atmospheres - 
Stars: Population II - globular clusters: general 



1. Introduction 

It is widely recognized that globular clusters (GCs) are 
cornerstones for the solution of a large variety of problems 
concerning the formation and evolution of galaxies. They 
are among the oldest objects formed in the Galaxy and, as 
a consequence, can be used as tracers of the early chemical 
evolutionary phases of the galactic environment. In princi- 
ple, as a first guess, one can consider the total metaLabun- 
dance (traditionally indicated by the ratio [Fe/H]EJ) of a 
globular cluster as representative of the original composi- 
tion of the gas from which it formed, and the variations in 
[Fe/H] among clusters as a fossil recording of the chemical 
enrichment history occurred in the Galaxy. It is nowadays 
well known, however, that this simple view is complicated 
by other phenomena: the chemical composition of the at- 
mospheres of stars can be altered by processes occurring 
during their evolution (see Smith 1987, and Kraft 1994 
for reviews). Even if the actual mechanisms are not com- 
pletely understood yet, it seems that surface abundances 
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^ We adopt the usual spectroscopic notation, i.e. [X] = 
log(X)star— log(X)0 for any abundance quantity X, and log 
e(X) — log (Nx/Nh) -I- 12.0 for absolute number density 
abundances. 
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of only the lightest element (C, N, O, Na, and Al) are af- 
fected, while Fc abundances arc unchanged. In a following 
paper we will discuss in detail this aspect. 

In the present paper we explore the possibility of build- 
ing a new mctallicity scale for galactic globular clusters us- 
ing only [Fe/H] values obtained from fine analysis of high 
dispersion spectra. Although globular cluster distances re- 
strict applicability of this technique to the brightest gi- 
ants, it provides direct accurate and quantitative deter- 
mination of the actual metallicity of stellar atmospheres. 

It is rather surprising that despite the advent of CCD 
detectors and of sophisticated and efficient spectrographs, 
the increasing number of measures of cluster metallici- 
ties has been (and is!) almost totally ignored in a variety 
of astrophysical problems involving this parameter. The 
most widely used metallicity scale for globular clusters is 
in fact the one obtained by Zinn & West (1984: hereinafter 
ZW) and Zinn (1985) from a calibration of integrated pa- 
rameters of globular clusters. The main advantage of us- 
ing integrated parameters is that they can be easily mea- 
sured even for the most distant objects in the Galaxy: 
homogeneous results can then be obtained for almost all 
known galactic clusters. However, integrated parameters 
are not directly related to metal abundances, and their 
use as abundance indices requires an accurate calibration 
in terms of the actual content of [Fc/H]. Reflecting un- 
certainties present at that epoch in abundances from high 
dispersion spectra, ZW attributed very low weight to di- 
rect abundance determinations for globular cluster stars 
when they constructed their metallicity scale for globular 
clusters. 

Since ZW work, a number of accurate high-resolution 
spectroscopic determinations of metal abundances for 
stars in globular clusters appeared in the literature. How- 
ever, only the work of Gratton and coworkers (Gratton et 
al 1986, Gratton 1987, Gratton and Ortolani 1989: G86, 
G87, G89 respectively; collectively G8689) was aimed at a 
systematic determination of abundances for a large num- 
ber of clusters (spectra for giants in 17 clusters were ac- 
tually analyzed). However, since completion of the Grat- 
ton and coworkers survey, there has been signiflcant pro- 
gresses both in the analysis techniques and on observing 
facilities. In fact the new Kurucz (1992, hereinafter K92) 
model atmospheres allow an homogcneoiis comparison be- 
tween solar and stellar abundances, a major drawback of 
former analysis of abundances for globular cluster stars 
(see e.g., Leep et al. 1987). Furthermore, improvements 
in high resolution spectrographs and detectors allow bet- 
ter spectra to be obtained for a larger number of stars: 
a major contribution has been done by observations with 
the Hamilton spectrograph at Lick by Kraft, Sneden and 
coworkers (Sneden et al. 1991; Kraft et al. 1992; Sneden 
et al. 1992; Kraft et al. 1993; Sneden et al. 1994; Kraft 
et al. 1995; hereinafter, SKPLl, SKPL2, SKPL3, SKPL4, 
SKPL5, SKPL6 and SKPL on the whole, for brevity). Re- 
gretfully, a homogeneous abundance scale based on high- 



resolution spectroscopic data does not exist yet, mainly 

due to inconsistencies in the model atmospheres and in 
the atomic parameters adopted in the various investiga- 
tions. Moreover, the need for such an improved scale is 
continuously growing, due to the high degree of accuracy 
required by a variety of problems, one for all: the long- 
debated calibration of the absolute magnitude of the hor- 
izontal branch in terms of the metal abundance. Minor 
variations in the adopted metallicities could result, ulti- 
mately, in non-trivial changes in globular clusters ages. 

Our goal is to exploit recent observational and the- 
oretical progresses to construct a new, reliable metallic- 
ity scale for globular clusters, completely based on high- 
quality data for red giant stars in 24 clusters. To this 
purpose, we obtained new data for a few stars in three 
southern clusters using the Long Camera mounted on the 
ESQ CASPEC spectrograph, and reanalyzed published 
equivalent widths (EWs), mainly from the Gratton and 
coworkers and Lick surveys, integrated by other sources of 
similar quality for clusters not included in those studies. 
These data were analyzed in a totally self-consistent way, 
allowing a modern calibration of the abundance indices 
considered by ZW. 

In Section 2 we present the new data, and in Section 
3 the data adopted from literature. In Sections 4 and 5 
we discuss respectively the atmospheric and atomic pa- 
rameters required for the abundance analysis. Our results 
will be exposed in Section 6, together with a discussion of 
the error sources. After a comparison with previous works 
(Section 7), we present our conclusions in Section 8. 



Table 1. Observed spectral regions 



Ac (A) Spectral interv. Indicators 



4200 


3800- 


-4600 


CH (G band) 


5100 


4700- 


-5480 


C2 


6300 


5950 


6700 


[0 I] red doublet 


8000 


7750- 


-8380 


ON red system 



2. Observational data 

We selected a sample of 10 red giant branch (RGB) stars 
in 3 GCs, representative of the typical range of metal 
abundance of these objects, i.e. 47 Tuc (high metallic- 
ity, [Fe/H]~ -0.8 dcx, 3 stars), NGC 6752 (intermediate 
metalhcity [Fc/H]~ -1.5, 4 stars) and NGC 6397 (low 
metaUicity, [Fe/H]~ -2.0, 3 stars). 

Due to obvious flux limitations, we restricted our ob- 
servations to the brightest globular cluster stars to ob- 
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Table 2. Program stars observed in 47 Tuc, NGC 6397 and 
NGC 6752 



Cluster Star V B-V {V - K)o 



47 TUC 


5422 


12.47 


1.40 


3.23 




5529 


11.87 


1.59 


3.86 




8406 


12.37 


1.43 


3.33 


NGC 6397 


C25 


12.22 


0.96 


2.29 




C211 


10.16 


1.46 


3.13 




C428 


11.50 


1.05 


2.47 


NGC 6752 


A31 


10.80 


1.60 


3.69 




A45 


11.57 


1.23 


3.00 




A61 


11.71 


1.13 


2.91 




C9 


12.37 


1.03 


2.62 



References - 47 Tuc: Frogel et al. (1981); NGC 6397 and NGC 
6752: Frogel et al. (1983). 



tain observational material good enough (i.e. adequately 
high resolution and S/N ratios) for fine abundance anal- 
ysis with reasonable exposure times. All observed stars 
were brighter than V=12.5; the faintest giant observed, 
star 5422 in 47 Tuc, has ^ -0.96 (adopting a true 
distance modulus of 13.31 and Ay — 0.12 from Djorgov- 
ski 1993). Moreover, stars were selected to have infrared 
photometry available, in particular in the K band, for an 
accurate determination of the effective temperature Tos- 

The observational material was acquired in two runs: 
October 1990 and June 1991. In both runs, echelle spectra 
of the program stars were obtained with the CASPEC 
spectrograph (in the Long Camera configuration) at the 
3.6m ESO telescope at La Silla, using a 31.6 lines/mm 
echelle grating. The slit width was adjusted in order to 
give resolution R ^ 30, 000. 

We tried to obtain a spectral coverage as large as pos- 
sible, to observe spectral features of different atomic and 
molecular species and to compare abundances of the same 
element as derived from different indicators. In Table [l] 
we list the main features (abundance indicators) of the 
spectral regions observed. Table || lists the literature pho- 
tometric measurements of the program stars, and Table 3 
the observed intervals for each star. 

Exposure times ranged from 15 to 70 minutes to re- 
duce cosmic ray contamination; we usually tried to ob- 
tain more than one spectrum for each object to eliminate 
the spurious events. A quartz lamp (for flat fielding) and 
a Thorium- Argon lamp (for wavelength calibration) were 
acquired after each program star exposure, with the tele- 
scope at the same position of the program star exposures. 
Besides, fast-rotating early- type stars were observed each 
night to remove telluric lines (see below). Bias frames have 
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been taken at the beginning of each night to account for 
readout noise. 



Table 3. Spectral intervals observed for the program stars 



Cluster 


Star 


Ac 
(A) 


^cxp 

(min.) 


S/N 


47 Tuc 


5422 


5100 


240 


145 






6300 


100 


111 






8000 


59 


45 




5529 


4200 


70 


60 






5100 


70 


100 






6300 


35 


56 






8000 


120 


100 




04U6 


OlUU 


z40 


1 on 
IzO 






6300 


90 


107 


NGC 6397 


C25 


6300 


80 


88 




C211 


4200 


40 


50 






6300 


15 


117 




C428 


6300 


41 


91 


NGC 6752 


A31 


6300 


30 


130 




A45 


6300 


45 


105 




A61 


6300 


50 


101 




C9 


6300 


90 


104 



2.1. Data analysis and equivalent widths 

The first steps in CCD reduction (bias subtraction, echelle 
order identification, scattered light subtraction, order ex- 
traction and wavelength calibration) were performed us- 
ing standard packages implemented in IRAFH environ- 
ment. Off-order scattered light was eliminated through 
bi-dimensional fitting along the dispersion and in the or- 
thogonal direction. The spectra were then wavelength cal- 
ibrated using a dispersion solution in two dimensions, de- 
rived from the Th-Ar lamps taken after each spectrum, 
and one-dimensional spectra were extracted using an op- 
timal extraction algorithm implemented in the package. 

The next steps of the analysis were then performed us- 
ing the ISA package (Gratton 1988), purposely developed 
to deal with high-resolution spectra. 

The blaze function was taken into account by divid- 
ing the spectra by dome fiat fields; the continuum was 
then traced on each individual spectrum for every ob- 
ject. Whenever we had multiple exposures of the same 

^ IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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star, spurious events and spikes due to cosmic rays were 
eliminated comparing different spectra. We then used the 
spectra of featureless, rapidly rotating, early-type stars 
for accurate removal of the telluric O2 features, affecting 
in particular the 6,300 A region. We first identified atmo- 
spheric features in the solar spectrum then we measured 
each line in the spectra of early-type stars acquired at dif- 
ferent airmasses z. From these measurements we derived a 
mean relationship between the EWs and z. For each star 
a synthetic spectrum of the atmospheric lines was then 
computed and convolved with the instrumental profile; O2 
features were finally cancelled out by dividing the spec- 
trum of each program star by the appropriate synthetic 
spectrum. This procedure allows to correct each star for 
the appropriate airmass; moreover it does not introduce 
additional noise in the object spectra. 

The resulting cleaned spectra were then added to im- 
prove the SjN ^ after correcting for the change in the (geo- 
centric) radial velocity of the star. Finally, a new contin- 
uum was traced. Figure ^ shows the tracing of a portion 
of our co-added and normalized CASPEC spectra for two 
program stars. 



■ ^ 




A7\f 




f 
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1 

84 
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— 1 



6290 6300 6310 



MA) 



Fig. 1. Tracing of a portion of the co-added, normalized spec- 
tra of 2 stars in the 6,300 A region. 

Final available spectra are listed in Table |^. We note 
that the following abundance analysis uses only EWs 
measured in the spectral regions centered at 5,100 and 
6,300 A. The main reason is that they are less affected 
by telluric bands and richer of stellar features (as com- 
pared to the 8,000 A region), and there is less concern 
related to line crowding and continuum level identifica- 



tion (as compared to the 4,200 A region); this last feature 
is of particular importance for the the spectra of stars in 
the metal-rich cluster 47 Tuc. 

Equivalent widhts of various element were then mea- 
sured in the two mentioned regions: in the following we 
will refer only to iron lines. Gaussian profiles were fitted 
to the observed profiles; when the number of clean lines 
was very low {e.g., for Fe II lines), we first derived a re- 
lationship between EW and central depth rc from un- 
blended lines and then we used it to add some new EW's 
for measured re's. The number of measured lines depends 
on S/N and on the star metallicity; generally, some 25-50 
Fe lines were measured for NGC 6397 stars, some 50-70 
for NGC 6752 stars and from 50 to 100 for 47 Tuc stars. 
In the following analysis, only lines with EW > 10 mA 
were used. Line parameters (see Section 5) and EWs are 
listed in Table |, both for Fe I and Fe II lines. 

We have 3 stars (namely, 47 Tuc-5529, NGC 6397- 
C211 and NGC 6752-A45) in common with another re- 
cent, high dispersion analysis by Norris and Da Costa 
(1995; hereinafter, NDC), at about the same resolution 
we used: this allows a comparison between the two sets 
of EW^, which is shown in Figure |[ The average differ- 
ence 1st i. EWcG96 - EWnbc = 6.1 ± 1.1 mA (ct=9.2 mA, 
76 linesa) . Note that here we regard this comparison as an 
external check of the accuracy of our EW^s measurements, 
but since we used data from NDC to enlarge the sample 
of analyzed stars (see below), the above comparison has 
to be regarded also as a self-consistency test on our total 
set of EWs. 

3. Literature data 

To increase our sample and to obtain a more statistically 
significant basis for our discussion, we also analyzed in a 
homogeneous way the globular cluster giants studied in 
two recent surveys: G8689 (41 red giants in 17 globular 
clusters) and SKPL (82 giants in 7 clusters). To these sets, 
we also added a full re-analysis of 3 giants in NGC 2298 
(from McWiUiam et al. 1992; McW92), 18 stars studied 
by Minniti et al. (1993; M93) in 8 clusters, and 8 stars 
in 2 clusters from NDC. Table || lists the whole sample of 
globular clusters red giants analyzed in the present paper, 
with the source for each subset of data. 

We simply adopted published EWs for all samples, 
since the observational material consists in high-resolution 
{R ~ 30, 000), high S/N spectra, fairly comparable to our 
own spectra. However, for the G8689 sample we had to 
make some corrections to the published EWs, before using 

We will indicate with the subscript CG96 the sample of our 
observed program stars and simply CG the results obtained 
using the whole sample (observations plus literature data) an- 
alyzed in the present paper 

^ Throughout this paper, the symbol a will indicate the stan- 
dard deviation of a single measurement, while the value after 
the symbol ± will refer to the standard deviation of the mean 
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Element A(A) E.P. log gf 47Tuc 47Tuc 47Tuc NGC6397 NGC6397 NGC6397 NGC6752 NGC6752 NGC6752 NGC6752 

5422 5529 8406 25 C211 C428 A31 A45 A61 C9 

EW EW EW EW EW EW EW EW EW EW 



Fe I 


5217.30 


3 


21 


-1 


07 


138 


6 


0.0 












































0.0 


Fe I 


5250.21 





12 


-4 


90 








0.0 


197 









































0.0 


Pe I 


5253.02 


2 


28 


-3 


81 








0.0 


85 


9 






































0.0 


Fe I 


5956.70 





86 


-4 


57 


166 


4 


0.0 


177 


9 


31 


3 


101 


3 


41 


7 








116 


8 


103 





0.0 


Pe I 


5976.79 


3 


94 


-1 


29 


92 


7 


0.0 


101 











40 


9 


24 


D 


70 


1 


55 





38 


7 


0.0 


Fe I 


5984.82 


4 


73 


-0 


28 








0.0 














35 





17 
















42 


3 


0.0 


Fe I 


6015.24 


2 


22 


-4 


66 


33 


5 


0.0 


49 


9 




















21 

















0.0 


Pe I 


6019.37 


3 


57 


-3 


21 


15 


4 


24.9 


15 


8 






































0.0 


Fe I 


6027.06 


4 


07 


-1 


18 


92 


6 


0.0 


94 


7 








29 


7 








69 

















0.0 


Pe I 


6034.04 


4 


31 


-2 


35 


14 





0.0 


27 


6 






































0.0 


Fe I 


6054.08 


4 


37 


-2 


22 


21 


5 


0.0 












































0.0 


Pe I 


6056.01 


4 


73 


-0 


45 


64 


6 


0.0 


81 


1 








18 


4 


11 


3 


52 


8 


43 


1 


37 


9 


0.0 


Fe I 


6078.50 


4 


79 


-0 


37 


76 


4 


0.0 


81 


4 






































43.9 


Pe I 


6079.01 


4 


65 


-0 


95 


54 


9 


0.0 


56 


5 




















22 


7 


32 


6 








0.0 


Fe I 


6082.72 


2 


22 


-3 


59 








0.0 














40 


4 














62 


7 








0.0 


Pe I 


6093.65 


4 


61 


-1 


29 


33 


3 


0.0 


41 


2 






































0.0 


Fe I 


6094.38 


4 


65 


-1 


52 








0.0 


17 








































0.0 


Fe I 


6096.67 


3 


98 


-1 


73 


53 


5 


0.0 


82 


3 




















41 


7 


29 


8 








0.0 


Pe I 


6098.25 


4 


56 


-1 


77 








0.0 


41 


9 






































0.0 


Fe I 


6120.26 





91 




94 


76 





95.4 


78 


1 




















64 


8 


32 


5 








0.0 


Pe I 


6151.62 


2 


18 


-3 


26 


124 


6 


133.5 


130 


2 








59 


6 


23 





125 


2 








57 


4 


57.7 


Fe I 


6157.73 


4 


07 


-1 


25 








0.0 


107 


6 








34 


5 








78 


8 








48 


5 


0.0 


Pe I 


6165.36 


4 


14 


-1 


46 


64 


5 


83.1 


71 


3 




















46 


4 


29 


7 








13.8 


Fe I 


6173.34 


2 


22 


-2 


85 


145 


9 


0.0 


147 


1 








82 


6 


46 


6 


157 


7 


94 


2 


91 


5 


87.2 


Pe I 


6187.40 


2 


83 


-4 


11 


33 


2 


0.0 












































0.0 


Fe I 


6187.99 


3 


94 


-1 


58 


69 


6 


0.0 


83 























47 


7 














23.5 


Pe I 


6213.44 


2 


22 


-2 


55 


157 


8 


0.0 


184 


6 


38 


8 


103 


4 


54 





167 




112 


7 


116 


3 


95.6 


Pe I 


6219.29 


2 


20 


-2 


42 


170 


2 


184.3 


180 





62 


8 


115 


8 


66 


5 


187 


9 


134 


2 


131 


1 


116.3 


Fe I 


6226.74 


3 


88 


-2 


02 








0.0 


83 


1 




















29 


6 














0.0 


Pe I 


6240.65 


2 


22 


-3 


23 








0.0 


137 









































0.0 


Fe I 


6280.62 





86 


-4 


35 


172 


1 


0.0 


203 


4 














56 


7 




















0.0 


Pe I 


6297.80 


2 


22 


-2 


73 


170 


9 


0.0 


180 


5 


38 


6 


99 


1 


50 


6 


152 


6 


113 


2 


112 


8 


90.6 


Fe I 


6301.51 


3 


65 


-0 


72 








0.0 








55 


7 














134 


6 














0.0 


Pe I 


6311.51 


2 


83 


-3 


14 








0.0 


























72 


D 


30 


D 


21 


7 


0.0 


Fe I 


6315.82 


4 


07 


-1 


65 


80 


4 


0.0 


90 


5 


























37 


6 








0.0 


Pe I 


6322.69 


2 


59 


-2 


39 


135 


3 


0.0 


149 


2 








79 


3 








141 


8 


104 


8 


93 


7 


72.6 


Fe I 


6353.84 





91 


-6 


44 


53 


6 


0.0 












































0.0 


Pe I 


6380.75 


4 


19 


-1 


32 


78 


7 


0.0 


90 


3 








16 


5 


15 


2 


52 


8 


40 





22 


6 


34.0 


Pe I 


6392.54 


2 


28 


-3 


95 


74 


3 


74.2 


69 


8 


























25 


D 


21 


7 


0.0 


Fe I 


6393.61 


2 


43 


-1 


43 


216 


1 


0.0 








87 


2 








99 


6 


222 


5 














0.0 


Pe I 
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them for a new analysis. In fact, in G87 the authors them- 
selves noted that their EWs were systematically higher 
than Mackle et al. (1975) EWs for a Boo (Arcturus). This 
is due to the lower resolution (R ~ 15, 000) used by G8689: 
unnoticed lines blend with the measured ones, resulting in 
an overall overestimate of the EWs. We derived the fol- 
lowing relation between G8689 EWs (EWq) and Mackle 
et al. (1975) EWs (EWm) for a Boo : 



EWm = 0.914 (±0.044)£;W^G - 14.3 (±17.6) niA (1) 

Strictly, this correction applies only to EWs measured on 
spectra having metallicities similar to that of Arcturus 
([Fe/H]~ —0.5); we expect smaller corrections for more 
metal-poor stars. To verify this point we considered the 
stars in common between G8689 and the present data 
(NGC 6397-C428 and NGC 6752-C9); in both cases there 
are 17 lines in common between the old and the new anal- 
yses. The mean difference between the EWs of our new 



analysis EWncw and those of G8689 for these two stars 
are: 

EWncw - EWg ^ -S.0± 1.8 ml (2) 

{a = 10.5 mA, 34 lines). Had we applied first eq. (|^) to 
correct the G8689's EWs, this difference would have been: 

EWnew - EWccorr = 15.3 ± 1.8 mA (3) 

(cr = 10.3 mA, 34 hues), where EW^corr are G8689 EWs 
after transformation using eq. (0). The overall correction 
is ~ 18 mA. We then estimate that while eq. (|l]) gives 
the appropriate corrections to G8689 EWs for metallic- 
ities similar to that of Arcturus ([Fe/H]^ —0.5), these 
corrections should reduce to ~ 1/6 of this amount for 
a metallicity in between those of NGC 6752 and NGC 
6397 ([Fe/H]- -1.7). We then included a metallicity- 
dependent correction, obtaining: 

EWc = EWg - (0.7 [Fe/H] -t^ 1.35) [(1 - a) EWq - b] (4) 
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_NCC6^97-Cail 



47TUC-5529 



NGC6752-A45 



Table 5. Complete list of globular clusters analyzed 



Fig. 2. Comparison of our EWs with those of NDC for the 3 
stars in common 





150 



I 1 I I 1 1 1 1 1 1 1 

_ B NGC6:597-C428 
. □ NGCf)75y C9 


1 1 1 






□ p 

'" 1 1 1 1 1 1 1 1 1 


(a) - 


_ pMlO-I-15 
_ X M71-I-31 
. B M71 -1-53 ^ 


II ,, 


- ^ /' 




□ M^^^ 

□ ,-' 

- . B/ ° 




!]-■□ □ 




" .ss 

/° , , , 1 , , , , 1 


(b) ■ 



GC 


Nr. Stars 


Reference 


NGC 104-47 Tuc 


3 


CG96 
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Fig. 3. (a) Comparison between EWs of G8689 sample after 
the correction and the EWs of the present analysis (sample 
CG96) for NGC 6752-C9 and NGC 6397-C428; (b) another 
test of self-consistency of our EWs, done by comparing the cor- 
rected EWs from G8689 sample with those from SKPL sample 
for 3 stars in common between the two sets. 



Eq. (|) was then applied to the old EWqs of G8689 
to bring them to the same system of our higher resolu- 
tion spectra. Figure^ (a) compares EWc with EWacw 
for the two stars in common between the G8689 and the 
CG96 sample; the average difference is EWc ~ EW^^^w = 
3 ± 2 mA (cr = 10 mA, 34 hues). The best fit fine 
(with constant term put to zero) to the data in Figure ^ 
(a) is: EWc = 1-025 ± QMbEW,,^^, with a scatter of 
11 mA. In the foUowing analisys, we used only lines with 
EWc >20 mA. 
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Figure ^ (b) shows another internal comparison, using 
G8689 (corrected) and SKPL EWs for 3 stars: the agree- 
ment is very good, the average difference being EWsKPh — 
-EWG8689 = -0.5 ±2.4 mA (cr=14 mA, 36 lines). We esti- 
mated internal errors in EWs of about ±9 mA in G8689 
set and ±5 mA in CG96, SKPL and the other sets ana- 
lyzed. We can then be confident that our results are based 
on an homogeneous and self-consistent dataset. 

4. Atmospheric parameters 

Model atmospheres appropriate for each star were ex- 
tracted from the grid of K92tl Values of effective temper- 
ature Tefi and surface gravity log g for giants observed in 
47 Tuc, NGC 6752 and NGC 6397 (sample CG96) were 
taken from Frogel et al. (1981, 1983). Frogel et al's Tea's 
are based on visual- near infrared colours (mainly V — K), 
transformed into temperatures based on the Cohen et al. 
(1978; hereinafter CFP) scale. \ogg values are determined 
from the position of the stars in the colour-magnitude di- 
agram (CMD). These values are considered as accurate as 
±100 K in Tes and ±0.3 dex in log g, including possible 
uncertainties in the effective temperature scale, as well in 
the adopted red giant mass, in the bolometric corrections 
and in the cluster distance modulus. 

Teff and log g values for stars in the SKPL sample are 
also generally based on CFP's scale, except for giants in 
M 92 and M 15; in these two cases the authors used Toff's 
from Carbon et al. (1982) and Trefzger et al. (1983). The 
reader should be warned that final values for Teff and log 
g, in the SKPL sample, are slightly different from the ac- 
tual values on the CFP's scale. This is because SKPL ad- 
justed the adopted photometric temperatures until a slope 
close to zero was achieved for the relationship between Fe I 
abundances from individual lines and excitation poten- 
tial. However, as discussed in their papers, the agreement 
of photometric and final spectroscopic Toff was from the 
beginning well within the uncertainties of the tempera- 
ture scale itself, so that neither the abundance ratios for 
individual stars nor the mean values for the clusters are 
significatively changed by adopting the tabulated values 
of atmospheric parameters from SKPL. 

In the G8689 sample, Teff and log g values were homo- 
geneously adopted from Frogel et al. (1979, 1981, 1983) 
and are also used in the present work. For other data sets 

^ The model atmospheres considered in this paper are those 
in Kurucz's CD-ROM 13. In these models, convection descrip- 
tion include an approximate consideration of possible over- 
shooting in the stellar atmospheres. This point is discussed at 
length in Castelli et al. (1996), who conclude that at present 
is not possible to establish whether or not this treatment of 
convection should be preferred to a more traditional approach 
without any overshooting. However, this uncertainty is of only 
minor concern for the metallicity scale established in this pa- 
per, insofar the approach used for the Sun and the globular 
cluster giants is the same. 



(McW92, NDC) we were able to use atmospheric param- 
eters from Frogel and coworkers. Only half of the 18 stars 
from M93 have Toff and logg values listed in Frogel et al.; 
for the missing stars, we used their original values, de- 
rived spectroscopically from the dependence of Toff on ex- 
citation potential and from the Fe ionization equilibrium. 
In fact, the parameters adopted by M93 are with very 
good approximation on the CFP's scale: the mean for the 
stars in common are: Toff (Frogel)-TefT(M93)=-29 ± 11 K 
(ct=30 K, 8 stars) and log g(Frogel)- log g(M93)=-0.10 ± 
0.06 dex {a— 0.17 dex, 7 stars). For further details, see the 
original papers of G8689 and SKPL. The starting input 
values for [Fe/H] were those from the original analyses. 

Metallicities were obtained varying the metal abun- 
dance [A/H] of the model until it was equal to the derived 
[Fe/H] value. For the microturbulent velocity Vt, the input 
values were changed until no trend in Fc abundance with 
the EW of Fe I lines was present Q. However, for giants in 
M 92 and M 15 only a few lines were available, so we took 
the values of vt from the relationship: 

vt = -0.322(±0.048) • log.g + 2.22(±0.31) (5) 

derived elsewhere (Carretta, Gratton & Sneden 1996; 
Gratton & Carretta 1996) for field stars. The adopted at- 
mospheric parameters are listed in Table 9 (also available 
in electronic form) for all the stars studied in the present 
paper. Metallicities were obtained varying the metal abun- 
dance [A/H] of the model until it was equal to the derived 
[Fe/H] value. 

5. Oscillator strengths, line selection and solar 
abundances 

Only lines clean from blends on very high resolution solar 
spectra were considered in the analysis; for Fe, the line list 
was extracted from Rutten & van der Zalm (1984), and 
Blackwell et al. (1980). The adopted values for the oscil- 
lator strengths gfs were determined following the same 
precepts of Clementini et al. (1995), who performed a 
high-dispersion study of the metal abundances of field RR 
Lyrae variables. Briefly, whenever possible laboratory gfs 
were considered: for Fe I lines they were taken from pa- 
pers of the Oxford group (for references, see Simmons & 
Blackwell, 1982) and Bard et al. (1991) and Bard & Koch 
1994), gfs of the Oxford group being corrected upward to 
account for the systematic difference with those of Bard 
et al. (0.03 dex; see Clementini et al. 1995). gfs for Fe II 
lines were taken from Heise & Kock (1990), Biemont et al. 
(1991) and Hannaford et al. (1992). For lines lacking ac- 
curate laboratory determinations, gfs were derived from 
an inverse solar analysis using the Holweger and Miiller 

® Expected line strengths were used when determining the 
microturbulent velocity, following the prescription of Magain 
(1984) 
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(1974: HM) model atmosphere and the Fe abundance de- 
rived from the other lines. The adopted gf values are also 
listed in Table ^ they can be integrated for other lines 
using gfs in the larger line list of Clementini et al. (1995, 
tables 3(a),(b)). 

In any differential analysis (i.e. a comparative anal- 
ysis in which the zero point of the [Fe/H] scale is set by 
the Sun) the assumed solar abundance is of paramount im- 
portance. The solar Fe abundance obtained with the set of 
gfs described above and the K92 solar model atmosphere 
(log e(Fe)=7.52), taken as the reference value for the Sun 
throughout the present work, agrees with the meteoritic 
value of Anders & Grevesse (1989: log e(Fe)=7.51). It is 
also very close to the value obtained using the HM model 
atmosphere (log e(Fe)=7.56: Castelh et al. 1996). A very 
similar value is obtained using Fe II lines when adopting 
the K92 model. The trend with excitation potential is also 
very small in the Sun. 

We wish to stress here the importance of the consis- 
tency between these various determinations of the solar Fe 
abundances: in fact, the whole scale of metal abundance 
previously determined for globular clusters was uncertain, 
due to the rather large difference (~ 0.15 dex) between 
the solar Fe abundances obtained using the HM and the 
Bell et al. (1976: hereinafter BEGN) model atmospheres 
generally adopted in the analysis of cool cluster giants (see 
e.g., Leep et al. 1987). In fact, it was not clear what solar 
Fe abundance to use: either the one determined using a 
model extracted from the same grid used for cluster gi- 
ants, or that obtained using the best solar model. Use of 
K92 atmospheres solves this problem, since we may now 
use models for giants extracted from the same grid which 
gives a solar Fe abundance in agreement with the best 
photospheric and meteoritic determinations. 

Finally, we note that this choice of lines and of gf 
values allows a direct comparison with the metallicity scale 
for the low-dispersion index AS in RR Lyrae stars derived 
by Clementini et al. (1995). 

6. Results and error estimates 

Our Fe abundances for both the original program stars 
and for all reanalyzed stars are listed in Table 9 (last 7 
columns )I3. Final [Fe/H] values adopted (last column) are 
those derived from neutral lines alone, since the number of 
Fe II lines with accurate EWs was often too small. When 
individual clusters are considered, our [Fe/H] values do 
not show any trend with Toff on the whole range 3800- 
4900 K (which approximatively corresponds to a range of 
about 2.5 in Mboi)- 

Our [Fe/H] values are systematically higher than those 
of the original analyses: the systematic difference is 0.12 ± 
0.01 dex (cr=0.08, 162 stars), as displayed also in Figure ^ 

All stars that were in more than one data sets have been 
independently re-analyzed, and then their [Fe/H] values aver- 
aged before computing the cluster's mean [Fe/H]. 
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Fig. 4. Comparison of new [Fe/H] values with those from the 
original analyses. Different symbols represent stars of different 
samples studied. 

This difference is mainly due to our use of K92 model at- 
mospheres for both solar and stellar analysis. In fact, in 
previous analyses {e.g., both G8689 and SKPL) the solar 
Fe abundances were obtained using the HM model atmo- 
spheres, which is ~ 150 K warmer than the BEGN mod- 
els in the line formation region. We notice here that rel- 
ative abundances {i.e. abundances obtained using model 
atmospheres from the same grid for both the Sun and the 
program stars) are almost insensitive to the grid adopted 
(differences are < 0.03 dex). In this respect, our analysis 
combines the advantages of both differential and absolute 
analyses, since our abundances are referred to the Sun, 
and we used a solar model extracted from the same grid 
of model atmospheres used for the program stars. 

Data concerning the Fe ionization equilibrium are 
shown in Table ^ which lists the mean differences between 
abundances derived from neutral and singly ionized lines 
of Fe. These values have been computed both for the to- 
tal sample and for the different sub-samples studied. From 
this Table we conclude that there is an excellent agreement 
between abundances derived from Fe I and Fe II, with no 
trend with Toff or [Fe/H]. The lack of any trend over the 
whole range of temperature is very important, since in the 
past some analysis {e.g., Pilachowski et al. 1983) claimed 
that a discrepancy was present between these 2 iron abun- 
dances in stars cooler than 4,300 K. The implication was 
that in the very upper red giant branch the usual Lo- 
cal Thermodynamic Equilibrium (LTE) assumption had 
to be released or, at least, carefully verified by statistical 
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Table 6. Mean differences Fe I - Fe II in globular cluster giants 



Sample Nr. Stars <Fe I - Fe II> a 



All 


147 


0.00 ± 0.01 


0.12 


CG96 


10 


-0.05 ± 0.04 


0.11 


G8689 


29 


-0.10 ± 0.03 


0.14 


SKPL 


82 


+0.02 ± 0.01 


0.08 


McW 


3 


+0.13 ± 0.04 


0.07 


M93 


15 


+0.11 ± 0.04 


0.15 


NDC 


8 


-0.01 ± 0.02 


0.06 



equilibrium computations. Our results, however, strongly 
confirms the recent study of Clementini et al. (1995; see 
also the footnote below) that pointed out that departures 
from LTE cannot significatively affect abundance analyses 
for stars cooler than RR Lyrae variables. 

For the SKPL sample it should be noted again that in 
their original papers both photometric gravities and Tgff 
values were purposedly changed to obtain a match of the 
two [Fe/H] abundances within 0.05 dex. 

Table shows the dependance of the derived abun- 
dances from uncertainties in the adopted atmospheric pa- 
rameters; this is obtained by re-iterating the analysis while 
varying each time only one of the parameters. To show 
how these sensitivities change with overall metal abun- 
dance, we repeated this exercise for both a metal-rich (star 
8406 in 47 Tuc) and a metal-poor giant (star A61 in NGC 
6752). 

Entities of variations are quoted in Table |^: these val- 
ues are larger than errors likely present in the adopted at- 
mospheric parameters. This will be shown in the following 
discussion, where we will try to provide reasonable eval- 
uations for the uncertainties in the adopted atmospheric 
parameters. To this purpose, we compared expected scat- 
ters in Fe abundances within individual clusters and differ- 
ences between abundances provided by neutral and singly 
ionized lines with observed values. Relevant data for this 
last parameter can be easily obtained from Table g3. For 
the reasons above mentioned, we omit from the following 
discussion the value from the SKPL sample and we con- 
centrate instead on the other mean differences, for which 
the standard deviation a = 0.11 + 0.15 represents the ran- 
dom errors contribution, and the error of the mean (0.01 
+ 0.04) the contribution due to systematic errors. 



Following the non-LTE analysis of Clementini et al. 1995, 
no significant departures from LTE are to be expected in RGB 
stars. The differences found in Fe abundances from Fe I and 
Fe II lines are thus likely to be interpreted as due to errors in 
the analysis and in the atmospheric parameters. 



6.1. Systematic errors 

The relevance of systematic errors is always difficult to 
reliably assess. We do not think there are serious concerns 
related to the adopted gf scale. On the other side, un- 
certainties due to the adopted model atmospheres may 
be large since various important aspects (like convection, 
molecular opacities, and horizontal inhomogeneities) are 
far from being adequately known. Large trends of Fe abun- 
dances with excitation have been obtained in the analysis 
of field metal-poor giants by Dalle Ore (1992), Dalle Ore 
et al. (1996), Gratton & Sneden (1994), and Gratton et 
al. (1996), when using both BEGN and K92 model at- 
mospheres. These trends suggest that currently available 
model atmospheres are not fully adequate for at least some 
metal-poor giants (see e.g. Castelli et al. 1996). While ab- 
solute abundances are quite sensitive to this source of er- 
rors, the comparison of relative abundances obtained with 
different model atmosphere grids (K92 and BEGN) sug- 
gests that our [Fe/H] values are not heavily affected. How- 
ever, our analysis should obviously be repeated once im- 
proved model atmospheres for metal-poor giants become 
available. 

We need to concern less about possible errors in the 
adopted temperature scale (in our case, the CFP one). In 
fact, were the Toff scale largely in error, we would expect 
a rather large difference between average abundances pro- 
vided by neutral and singly ionized Fe lines. The values 
listed in Column 2 of Table indicate that a systematic 
error of 100 K in the adopted Teg's would translate into 
a systematic difference of 0.2 dex between abundances of 
Fe I and Fe II. Since the observed difference ranges from 
0.02 dex to 0.13 dex (depending on the considered sam- 
ple), we conclude that the Toff scale cannot be systemati- 
cally incorrect by more than 50 K. 

6.2. Internal errors 

Internal errors may be determined from a comparison with 
the observed scatter in our abundance determinations (of 
individual lines and of individual stars in each cluster). We 
will consider only errors in the EWs and in the adopted 
atmospheric parameters, while we regard internal errors 
in the adopted gfs as negligible. 

6.2.1. Equivalent widths 

The scatter of abundances from individual (Fe I) lines is 
0.13, 0.11, 0.15, 0.15, 0.14 and 0.12 dex for the CG96, 
SKPL, G8689, McW92, M93 and NDC samples respec- 
tively. These values for the scatter can be ascribed to er- 
rors in the EWs of a few mA (see Section 3), and yield 
mean internal errors of 0.03 and 0.06 dex for Fe I and Fe II 
respectively. These internal errors can be added quadrati- 
cally and give a prediction of about 0.07 dex for the scatter 
in the differences between abundances derived from Fe I 
and Fe II lines. Since the observed scatter ranges from 



10 E. Carretta, R. G. Gratton: Abundances for Globular Cluster Giants: 

Table 7. Dependence of the derived abundances on atmospheric parameters 



Element AT^a A log 3 A[A/H] Avt Atot Atot Amod 

+ 100K +0.5 dex +0.2 dex +0.5 km s"^ random syst. BEGN-K92 



Star 8406 in 47 Tuc ([A/H] = -0.61) 
[Fe/H]I 0.015 0.118 0.044 -0.232 0.100 0.080 -0.124 

[Fe/H]II -0.185 0.318 0.091 -0.094 0.136 0.183 -0.165 



Star A61 in NGC 6752 ([A/H] = -1.61) 
[Fe/H]I 0.149 0.017 -0.015 -0.094 0.082 0.079 -0.039 

[Fe/H]II -0.077 0.235 0.058 -0.053 0.081 0.117 -0.017 



a — 0.11 to a ~ 0.15 (depending on the adopted sample), 
additional sources of errors are clearly present, probably 
related to the adopted values for the atmospheric param- 
eters (see Table ^ and discussion below). 



\ogg and in Toft are not tiedB. But, as matter of fact, 
temperature and gravity are not completely independent, 
since to derive \ogg from the position of the star in the 
CMD we have to use the relationship L — AirR^aTcs'^, i.e.: 



6. 3. Temperatures 

CFP V — K colours have errors of ^ 0.05 mag, which cor- 
responds to 35-40 K using their calibration. This is the 
internal error of T^s's for stars within a cluster. When 
comparing stars in different clusters, the effects of errors 
in the interstellar reddening should also be considered. 
Comparing various estimates for the same cluster, we es- 
timate an uncertainty of 0.02 + 0.03 mag in E{B — V), 
and 2.7 times larger in E{V — K). Hence, there is an ad- 
ditional systematic error of ~ 0.05 mag in the {V — K)o 
colour 35 — 40 K) systematic for all stars in a clus- 
ter (but random from cluster to cluster) due to errors in 
the reddening. If we add these two uncertainties quadrat- 
ically, we estimate that the adopted Toff's have internal 
errors of ~ 50 K. The same figures approximately hold for 
the B — V colour, which is a less accurate temperature 
indicator (see e.g., Gratton et al. 1996), but at the same 
time is measured with a precision better by a factor of 5 
than the V — K for bright globular cluster giants. 

Table |^ suggests that most of the residual scatter in 
the differences between Fe I and Fe II abundances may be 
attributed to random errors in the adopted Toff values. 



6.4. Gravities 

The adopted gravities were deduced from the location of 
the stars in the CMD. Since they were not deduced from 
the ionization equilibrium, one could think that errors in 



log g/gQ = 4 log T,s/Tq - log T/Lq + log M/Mq (6) 

To estimate the order of magnitude of the errors affecting 
gravity, consider the following: 

— the luminosity log L can be wrong either if the clus- 
ter distance moduli or the bolometric corrections EC 
are wrong. The distance moduli were derived by CFP 
from a compilation of literature data; however uncer- 
tainties are not larger than 0.2 mag, i.e. 0.08 dex in 
log L. The error in the ECs depends on the model at- 
mospheres. Toff's, and metal abundances. It cannot be 
larger than 0.2-0.3 mag, though, else the luminosity 
at the tip of the RGB would disagree with that pre- 
dicted by stellar evolution models. This contribution 
then translates into another 0.08 dex. If we sum these 
two contributions quadratically, we estimate that the 
total uncertainty in adopted luminosities is ~ 0.11 dex 
in logT. 

— the contribution of errors in the mass is very small, 
since this is fixed by the age of the globular clusters, 
and it cannot be wrong by more than 10%. The uncer- 
tainty in logM may then be taken as 0.04 dex. 

— as mentioned above, uncertainties in Toff's are of ~ 
50 K, i.e. ~ 0.005 dex in log Toff. 

From these considerations, we estimate that the adopted 
gravities have internal errors of ^ 0.15 dex. 

In column 3 of Table ^ we investigate the effects of a 
variation of 0.5 dex in the surface gravity; on the basis 

' As a generic estimate, an error of 100 K in T^g translates 
into a 0.3 dex error in logg, when deriving spectroscopically 
the gravities 
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of the previous discussion, the contribution from this col- 
umn should be then divided by at least a factor of 3. It 
is interesting to note that a larger error of Alogg — 0.25 
would explain the whole residual 0.11 dex in the random 
error. This is not the case, though, since there is surely a 
contribution from errors in Tefj: this further confirms that 
A log 3 = 0.25 is an overestimate, and the assumed value 
of 0.15 dex is reliable. 

6.5. Metallicities 

For each star analyzed we have also random errors in the 
estimate of [A/H] due to errors in T^g, in gravity (of little 
entity) and in the measured EWs. This kind of errors 
can be evaluated from independent analyses of the same 
star. To this purpose, we can compare the results obtained 
for stars in the same cluster, since they are thought to 
share the same overall metallicity: the r.m.s deviation from 
the mean will give an idea of the uncertainties due to 
random factors. The quadratic average is 0.06 dex and so 
they contribute very little to the observed difference in the 
abundances from Fe I and Fe II (less than 0.025 dex, from 
Table 0). 

6.6. Microturhulent velocity 

The internal error in the vt is usually estimated from the 
comparison of empirical and theoretical curve-of-growth; 
it is typically not larger than 0.2 km s""'^ for the giants 
analyzed, since the microturhulent velocity is derived us- 
ing Fe I lines both on the linear and saturation part of 
the curve-of-growth. As above, an independent test of the 
random errors comes from the comparison between the 
values obtained for the same star independently analyzed. 
We obtained Awt=0.17 km s^^ for the star C428 in CG96 
and in the G8689 sample; it confirms that the microtur- 
hulent velocity has an error smaller than 0.2 km s^^. 

6. 7. Discussion of errors 

To conclude, we have to consider two kinds of errors: first, 
the internal, random errors, that affect the comparison 
from star to star, and second, the systematic errors, that 
give an idea of the reliability of our metallicity scale, of 
the temperature scale adopted, etc. For the random er- 
rors, we have seen that reasonable estimates are 50 K in 
Teff, 0.15 dex in logg, 0.06 dex in [A/H] and 0.2 km s^^ 
in Vt ; these errors will affect the scatter of our data. As to 
systematic errors, we have only the indication given by the 
difference in the abundances from neutral and singly ion- 
ized Fe lines; from the previous discussion, we conclude 
that these errors are of the same order of magnitude of 
random ones. 

Columns 6 and 7 of Table list the uncertainties in 
the [Fe/H] ratios derived from the quadratic sum of the 
contributions from random and systematic errors, respec- 
tively. We remark that the changes in the parameters used 



to construct these columns are not those indicated in the 
Table, but the more realistic estimates obtained from the 
above discussion. From Table |^ we can estimate that the 
total uncertainty in our Fe I abundances (from which we 
derive the clusters metallicity) is about 0.11 dex for the 
most metal-poor stars, increasing to about 0.13 dex for 
the most metal-rich stars. 

7. The metallicity scale of the globular clusters 

Mean values of the metallicities derived for the 21 clusters 
analyzed are listed in Table ^. 

The internal uncertainty in [Fe/H] abundances 
(a/N^^^, where N is the number of stars studied in each 
cluster) is very small: on average, 0.06 dex, which can be 
interpreted also as the mean precision of the cluster rank- 
ing on our new metallicity scale. For comparison, in the 
same Table, we also give the original [Fe/H] ratios ob- 
tained in previous analyses. In the last column the metal 
abundances from the compilation of Zinn and West (1984) 
are listed, superseded and integrated for a few clusters by 
the new measurements of Armandroff and Zinn (1988); 
this scale will be indicated as a whole, hereinafter, as ZW. 

7.1. Comparison with the ZW scale 

The 24 clusters of Table ^ can now be regarded as stan- 
dard reference clusters to calibrate individual metal abun- 
dance indicators with metallicities directly derived from 
high-dispersion spectroscopic analysis. We feel confident 
that our list covers fairly well the whole range in metal- 
licity of globular clusters, going from typical metal-rich 
clusters, as 47 Tuc, M 71 and NGC 6352, to the classical 
mctal-poor templates (M 92, M 15, M 68). The sample 
of intermediate metallicity clusters is also very well repre- 
sented among our calibrators. One of our main purposes 
is to revise and refine the calibration of the ZW ranking 
system, which covers almost all known globular clusters. 
The main advantage of the ZW scale is that their system 
is applicable even to the most distant objects, being based 
on the integrated parameter Q39 and/or on low-dispersion 
spectroscopy of the infrared Ca II triplet. On the other 
side, any integrated index is not, by definition, a function 
of a single element in a globular cluster. In particular, the 
major contribution to line blanketing in the spectral range 
covered by the Qsg index is due to the H and K lines of 
Ca II, with other significant fractions due to the A 3883 
CN band and some Fe blends. Hence, reliability of ZW 
metallicities ultimately rests on the coupling between Ca, 
C, N, and Fe abundances. It is outside the purposes of 
this paper to proceed further on this point; we only wish 
to recall here that the strength of CN-bands is known to 
vary from star to star (the so-called CN-signature) , having 
a bimodal distribution in most (but not all!) clusters (see 
Kraft 1994 for a recent review). 
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Table 8. Mean metallicities for globular clusters compared to literature data 



NGC Messier Stars Meanir.m.s. a G8689 SKPL Others ZW 



104 


47 Tuc 


5 


-0.70±0.03 


0.07 


-0.82 






-0.71 


288 




2 


-1.07±0.03 


0.04 


-1.31 






-1.40 


362 




1 


-1.15 




-1.18 






-1.27 


1904 


M 79 


2 


-1.37±0.00 


0.01 


-1.42 






-1.68 


2298 




3 


-1.74±0.04 


0.06 






-1.91" 


-1.85 


3201 




3 


-1.23±0.05 


0.09 


-1.34 






-1.61 


4590 


M 68 


3 


-1.99±0.06 


0.10 


-1.92 




-2.17*" 


-2.09 


4833 




3 


-1.58±0.01 


0.01 


-1.74 




-1.71'' 


-1.86 


5272 


M 3 


10 


-1.34±0.02 


0.06 




-1.46 




-1.66 


5897 




2 


-1.59±0.03 


0.05 


-1.84 






-1.68 


5904 


M 5 


16 


-1.11±0.03 


0.11 


-1.42 


-1.17 




-1.40 


6121 


M 4 


3 


-1.19±0.03 


0.06 


-1.32 






-1.33 


6144 




1 


-1.49 








-i.sg*" 


-1.75 


6205 


M 13 


23 


-1.39±0.01 


0.07 




-1.49 




-1.65 


6254 


M 10 


15 


-1.41±0.02 


0.10 


-1.42 


-1.52 




-1.60 


6341 


M 92 


9 


-2.16±0.02 


0.05 




-2.25 




-2.24 


6352 




3 


-0.64±0.06 


0.11 


-0.79 






-0.51 


6362 




2 


-0.96±0.00 


0.01 


-1.04 






-1.08 


6397 




10 


-1.82±0.04 


0.10 


-1.88 






-1.91 


6656 




3 


-1.48±0.03 


0.06 


-1.56 






-1.75 


6752 




12 


-1.42±0.02 


0.08 


-1.53 






-1.54 


6838 


M 71 


13 


-0.70±0.03 


0.09 


-0.81 


-0.79 




-0.58 


7078 


M 15 


4 


-2.12±0.01 


0.01 




-2.30 


-2.23*" 


-2.17 


7099 


M 30 


2 


-1.91±0.00 


0.00 






-2.11'' 


-2.13 



References: a = McW92, b = M93 



Moreover (see e.g., Clementini et al. 1995, section 
5.1.1) the [Ca/Fe] ratio does not scale with Fe on the 
whole range of metaUicities, being lower in metal-rich than 
in metal-poor Population II stars. Furthermore, a serious 
caveat has been advanced on the claimed independence 
of the Qag index from the horizontal branch morphology 
(see e.g., Smith 1984). To overcome this kind of problems, 
the most straightforward way to correct the ZW scale con- 
sists in working directly on the final metallicities, since the 
original compilation of Zinn and West (1984) was obtained 
averaging a number of [Fe/H] values derived from differ- 
ent indicators {e.g., {B — V)Q^g, Ayi.4, AS") and calibrated 
against Qsg- 

In Figure |^ we then compare our high-dispersion 
[Fe/H] values with the ZW values for the 24 calibrating 
clusters. The error bars (1 a) are from Zinn and West 
(1984: Table 5) and from our Table |^. As it is evident 
from this figure, the ZW scale is far from linear, deviat- 
ing both in the low and in the high metallicity regimes, 
when compared with [Fe/H] from our direct analysis. In 
the metal-rich region ([Fe/H]> —1) ZW's metallicities 
are on average 0.08 dex too high for the 3 clusters 47 
Tuc, M 71 and NGC 6352, with the last two objects 
being responsible for most of the discrepancy (0.12 and 



0.13 dex, respectively). For -1 <[Fe/H]< -1.9 the [Fe/H] 
values of ZW are definitively too low by a mean value 
of 0.23 dex (cr=0.09 dex, 16 clusters). Finally, in the very 
low-metallicity tail, ZW's values are on average 0.11 dex 
higher than ours (ct=0.06 dex, 6 clusters). 

The non-linear behaviour has been confirmed by a 
t— test on the significance of the quadratic term in the 
relation between ZW and ours [Fe/H]'s. To bring ZW's 
[Fe/H] ratios on a metallicity scale fully based only on high 
dispersion spectroscopy (HDS) we then derived a correc- 
tion given by a quadratic relation. This procedure auto- 
matically takes into account also the different zero point 
between the two scales, since the ZW scale was ultimately 
based on the Cohen (1983) scale, which, as other past 
analysis, adopts the traditional old solar Fe abundance 
log e=7.67. The resulting function we derive for this cor- 
rection is: 

[Fe/H]cG = -0.618(±0.083) 

-0.097 (±0.189) [Fe/H]zw 

-0.352 (±0.067) [Fe/H] zw^ (7) 

with the correlation coefficient r = 0.982 and f7=0.08 for 
24 clusters. This relationship is highly significant, from a 
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Fig. 5. Mean metallicities for the 24 clusters from the present 
work compared with metallicities on the Zinn and West scale 
(1984). 

statistical point of view, and can be applied to ZW metal- 
licities in the range —2.24 <[Fe/H]zvv < —0.51, defined by 
the lowest and highest values of [Fe/H]zw among the clus- 
ters used for the calibration. The quadratic regression line 
is shown as a heavy line in Figure overimposed in the 
same Figure is also the result of a linear fit, which takes 
into account the errors. As one can see, even considering 
3(7 error bars, it it very difficult to represent the data on 
a linear scale, in particular at the lower metallicity edge. 

Once the correction is applied the non-linearity of the 
ZW scale obviously disappears. However, a certain amount 
of scatter is still present in the intermediate-metallicity 
regime; we believe that it could be attributed to a resid- 
ual effect, not well removed by our calibration, of the sec- 
ond parameter. This last is in fact likely to affect ZW's 
metallicities more severely in this regime, in which the in- 
tegrated colours of clusters of different HB morphological 
type can be sensibly misinterpreted in terms of [Fe/H]. 

The next logical step would be now to calibrate other 
empirical metallicity indicators, i.e. repeat the original 
work of Zinn and West (1984) but using now our direct 
[Fe/H] values from HDS as a calibrating sequence. The 
most interesting and accessible parameters are the photo- 
metric ones [e.g., (B-V)o,g, AV1.4, etc.): they are widely 
used since it is easy enough to measure them from the re- 
cent and accurate CCD-based colour-magnitude diagrams 
(CMD). However, it would be preferable to have a dataset 
of homogeneity and accuracy comparable with the preci- 
sion of our metallicities, instead of relying on compilation 



from different sources. Since such an effort is presently in 
progress on a set of CMDs analyzed in a self-consistent 
way, we postpone to a forthcoming paper this kind of cal- 
ibration. However, an immediate and meaningful compar- 
ison can be made with the metallicity scale for globular 
clusters derived from RR Lyraes, since we can compare re- 
sults obtained for two different stellar populations, RGB 
stars and HB stars, independently checking the validity of 
both scales. 

7.2. Comparison with the metallicity scale of RR Lyrae 
stars 

The most recent calibration of [Fe/H] in terms of the Pre- 
ston's (1959) index AS is the one defined by Clementini 
et al. (1995), who found 

[Fe/H] = -0.194(±0.011)AS - 0.08(±0.18) (8) 

This relation was derived using RR Lyraes both in the 
field and in globular clusters. However, while metallicity 
values for field RR Lyraes were directly derived from high- 
resolution spectra or from the re- analysis of literature data 
(for a total of 28 RR Lyraes), cluster metallicities were 
taken at face value from the literature, even if a zero point 
was admittedly noted while using data from different sam- 
ples. We have many clusters in common with the study of 
Clementini et al. (1995) and have then derived again the 
[Fe/H] vs AS relation. Figure^ shows the result of our 
re- analysis. 

We obtained AS" values for 15 of our calibrating clus- 
ters from the metallicities of Costar & Smith (1988), in- 
verting the Butler's (1975) relation they used. Our values 
are not completely identical to those used by Clementini 
et al. (1995); the main differences are that a) we excluded 
47 Tuc, since its value for AS* is based on a single star, 
possibly not member of the cluster (Tucholke 1992) and 
b) for NGC 288 we assumed a value of AS* — 5.9, since 
the mean value 7.2 cited by Costar & Smith was obtained 
including spectra taken at phases near maximum light. 

Regression lines were then obtained by least-squares 
fits (we averaged values obtained exchanging the indepen- 
dent and dependent variables): 

— If we consider only the 15 GCs we obtain: 

[Fe/H] = -0.165 (± 0.019) AS* - 0.142 (±0.033) (9) 

with a o- = 0.130, and r = 0.947 (dotted fine in Fig- 
ure H). 

— It we consider both the 15 GCs and the 28 field RR 
Lyraes we get: 

[Fe/H] = -0.187 (±0.011) AS* - 0.088 (±0.041) (10) 

with a (T = 0.269, and r = 0.954 (solid line in Figure H). 

Also shown in Figure ^ is the calibration obtained by 
Clementini et al. (1995), using only 28 field variables (their 
equation 6): [Fe/H]= -0.204(±0.012)AS -0.102(±0.036), 
CT=0.190. 
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Fig. 6. Calibration of the AS index with our new analysis and 
with the data of field RR Lyrae variables from Clementini et 
al. (1995) 

The first striking evidence both from Figure ^ and the 
above equations is that the sequence of the globular clus- 
ter points seems to be much better defined, with a smaller 
scatter than the distribution of field RR Lyraes. The scat- 
ter in equation (6) of Clementini et al. is 0.19 dex, to be 
compared with the value of 0.13 dex obtained using only 
the new values for the clusters. We stress the fact that 
both the solar Fe abundance and the source for the oscilla- 
tor strengths are in common between the present analysis 
and that of Clementini et al. (1995); moreover, the pro- 
cedure followed in the abundance analysis is virtually the 
same. This may be evidence in favour of a larger intrinsic 
scatter in field than in cluster variables, or it may just 
reflect a smaller error in the values of AS" for cluster RR 
Lyraes. However the last explanation seems a little un- 
palatable, since determinations of AS* values are usually 
more accurate for nearby field stars. 

The second feature shown in Figure ^ is a rather clear 
separation between the relations for cluster and field RR 
Lyraes in the low metallicity region; this is the likely ex- 
planation for the increase in the scatter when the cali- 
bration [Fe/H]-AS is made using both cluster and field 
variables. The same behaviour was evident also in Figure 
14c of Clementini et al. (1995), but here it is even clearer, 
given the high degree of homogeneity in our data. Why 
this is so, we are not sure, apart from a suggestion of non- 
linearity in the AS'-[Fe/H] relation theoretically predicted 
(Manduca 1981) and discussed in Clementini et al. (1995). 
Apart from this, there seems to be a good agreement be- 



tween both scales; if we use our new calibration (equation 
12) to derive [Fe/H] ratios, the differences [Fe/H]cG — 
[Fe/H]AS are on average 0.12 ± 0.03 (cr=0.10, for 16 clus- 
ters). 

8. Discussion and conclusions 

We think that the effort to get more reliable and accurate 
metallicities is truly worthwhile: a new, homogeneous scale 
of [Fe/H] for globular clusters is really needed, since up to 
now the numerous but still heterogeneous estimates from 
HDS have been systematically ignored in many problems 
of stellar evolution. 




3,85 3.8 3.75 3.85 3.8 3.75 

l°g T.„ log T,„ 



Fig. 7. Observed distributions of variables in M 3 (filled 
squares) and M 15 (open squares) in the log P'- log Tcff plane 
(see text). 

As an illustrative example we will explore the effects 
of our new [Fe/H] values on the long-standing problem 
of the Oosterhoff (1944) effect among globular clusters, 
which belong to 2 groups on the basis of the mean period of 
their type ah RR Lyraes; this division reflects a separation 
in metallicity of the clusters (Arp 1955). The statement of 
this problem, its history and references are fully addressed 
in a recent paper by Sandage (1993a). Briefly, the concept 
is to use the pulsation equation for RR Lyrae stars (P 
— f (Tcff ,L,M), where P is the fundamental period of the 
pulsation, L the luminosity and M the mass of the star) 
to derive a calibration of the absolute luminosity of RR 
Lyraes in terms of the metal content, [Fe/H], in turn a 
cornerstone to ultimately get the ages of globular clusters. 
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This is possible if the variations of the parameters involved 
(P,Toff, L, M) with [Fe/H] are known. 

The critical and more debated point is to determine 
how AlogP, the shift existing in the log P-log Toff plane 
between the distribution of variables in Ool and OoII clus- 
ters, varies with the metallicity; the slope, in particular, 
is still much controversial. While Sandage (1993a) finds 
a value of about 0.12, theoretical models seem to predict 
a much lower value, around 0.05 (in the sense of longer 
periods for metal-poorer clusters). Our approach is based 
both on our new metallicities and new Tos's for variables 
in M 3 (Ool) and M IL. (OoII), i.e. the template-pair for 
the Oosterhoff effect. 

(a) From the present work we adopt [Fe/H]= —1.34 for 
M 3 and [Fe/H]= —2.12 for M 15. These values compare 
well to those derived by Sneden et al. (1991) and Kraft 
et al. (1992, [Fe/H]= -1.48 and -2.30 for M 3 and M 15 
respectively), allowing for their use of the old Bell et al. 
(1976) model atmospheres and of a different set of gf. 

If one compares these new values with the classical ZW 
ones (—1.66 for M 3 and —2.15 for M 15) , it it immediately 
evident that whatever the period shift between M 3 and M 
15 may be, it has to be "diluted" on a larger range of 
metallicity than before. This in turn affects the slope of 
the relationship log P-[Fe/H]. 

b) We adopt here the new temperature scale derived 
by Gratton et al. (1996). Briefly, they first derived empir- 
ical colour- Toff relations for population I stars, based on 
Teff's from the Infrared Flux Method (IRFM, BlackweU 
and Lynas Gray 1994; Bell and Gustafsson 1989) and in- 
terferometric diameters (Di Benedetto and Rabbia 1987); 
stars of luminosity class III and V were considered sep- 
arately. Toff's for stars of different log g and [Fe/H] are 
then obtained applying systematic corrections to the the- 
oretical K92 Toff's, to tie them to the empirical calibra- 
tions and to the K92 colours, in order to take into account 
the real metallicities (different from solar). The effect of 
gravity was taken into account by interpolating between 
typical values for dwarfs and giants. The underlying as- 
sumption is that K92 models (the same consistently used 
in the abundance analysis) are well able to reproduce the 
run of colours with the overall metal content [A/H], but 
have to be corrected in function of the surface gravity (a 
constant mean value of log g — 2.75 was assumed for all 
variables). 

The adopted scale is very similar to the one defined 
by Clementini et al. (1995) to study field RR Lyrae stars; 
they also showed that as far as colours (and abundances) 
are concerned, K92 models are well suited to reproduce 
the atmospheres of RR Lyrae variables.EJ 

Our conclusions would not change had we used M 68 
([Fe/H]= —1.99 on our scale) instead of M 15 as OoII tem- 
plate. While M 68 has a small and well determined reddening 
{E{B — V) =0.03), M 15 has a better populated instability strip. 

The reader should be aware of the small inconsistency due 
to our use of the CFP scale for globular cluster (cooler) giants 



The photometry for the variables of M 15 is taken from 
the high quality work of Bingham et al. (1984). For M 
3, we used colours from Sandage (1990). The choice of 
the mean colour which better represents the one the vari- 
able should have were it not pulsating is not simple (see 
e.g., Fernley 1994). However, for sake of comparison with 
Sandage's previous works on M 3, in the present study we 
use {< B > — < V >)corr, the corrected colours taken 
from the quoted sources; they also include an empirical 
correction AC (A) depending on the light-curve amplitude 
(Sandage 1990). 

We adopt E{B - V")=0.00 for M 3, while for M 15 we 
assume 0.10 (Zinn 1980). 

Figure shows the observed distributions of variables 
in M 3 and M 15 in the classical log P'~log Toff plane (log 
P' is the reduced and fundamentalized period; for c-type 
RR Lyraes the period is fundamentalized by adding 0.127 
to log P). The lines are best-fits drawn by eye through 
the data, since the computed least squares linear regres- 
sion result in non-parallel lines, due to the scatter among 
the data, especially in those of M 3. This method is good 
enough in this case, since we are only interested in pre- 
senting and emphasizing another source of uncertainty af- 
fecting this method. 

In panel (a) Toff's are derived from the old Kurucz 
(1979) model atmospheres, following the (unpublished) 
calibration of Buser, and with E{B - V)=0.00 and 0.10 
for M 3 and M 15 respectively. Metallicities are on the 
Zinn and West's scale: this panel should then reproduce 
closely enough Sandage's results for the Oosterhoff effect 
in the pair M 3-M 15. In panel (b) we used the new metal- 
licities from the present work and Toff's derived from K92 
models. 

In order to evaluate the period shift between the two 
distributions, we followed the prescriptions of Sandage 
(1993b), reading the values of log P' not at constant log 
Toff, but at lower temperatures for lower metallicities, fol- 
lowing his equation (5). Our measurements then give the 
period shift AlogP (in the sense M 15-M 3) = -0.103 and 
—0.076 for case (a) and (b), respectively. We note that 
only the first slope is similar to the one derived by Sandage 
(1993a: AlogP/A [Fe/H]= -0.122). Simply using the new 
metallicities and effective temperatures, although derived 
from the same photometric data, significatively decreases 
the size of the effect. 

Taking at face value these figures, and using the pa- 
rameters temperature and mass (along with their varia- 
tions with metallicity as given by Sandage 1993b) in the 
pulsation equation of van Albada and Baker (1973), we 
obtain for the slope of the relationship Mboi vs [Fe/H] the 
values 0.300 and 0.221 for case (a) and (b) respectively. As 



and the new temperature scale for RR Lyraes; however, as 
discussed in Gratton et al. (1996), differences between the two 
scales are well within the intrinsic uncertainties of both scales 
and do not alter significatively our conclusions. 
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one can see, case (b) is very similar to that derived from 
the Baadc-Wcssclink analyses (0.25, as quoted for instance 
in Sandage 1993b), whose results were considered up to 
now in serious disagreement with the calibration from the 
Oosterhoff effect. 

This straightforward exercise points out that simply 
using new, modern, and self-consistent temperatures and 
metallicitics, the size of the Oosterhoff effect in the tem- 
plate pair M 3-M 15 is somewhat reduced, being more 
consistent with theoretical models of Zero Age Horizon- 
tal Branch (see Sandage 1993b for references). We then 
conclude that our effort in obtaining these new, improved 
values for metallicities can pay off; we will proceed to a 
new calibration of photometric / low resolution indices, 
including all known globular clusters, as soon as a homo- 
geneous dataset will be available. 
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Table 9. Adopted atmospheric parameters and results of abundance analysis for globular cluster giants 
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